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The minimal supersymmetric standard model (MSSM) is extended by including an additional 
pair of constrained Higgs doublet superfields through which the electroweak symmetry is nonlinearly 
realized. The superpotential couplings of this constrained Higgs doublet pair to the MSSM Higgs 
doublet pair catalyze the later nontrivial vacuum expectation values. The mass spectrum of the 
Higgs scalars and Higgsino-gaugino sector is presented for several choices of supersymmetry (SUSY) 
breaking and Higgs superpotential mass parameters. The effect of the the additional constrained 
fields and multiple vacuum expectation values on the lowest mass neutral Higgs scalar production 
and decay is presented. 



As experimental results continue to squeeze the 
MSSM, particularly through the increasing lower 
bound on the Higgs boson mass, it proves worth- 
while to focus increasing attention on SUSY models 
beyond the MSSM. Indeed, there are a plethora of 
such models already introduced, with the NMSSM[T] 
garnering the most attention. Here we focus on an- 
other model (actually a class of models) |2j character- 
ized by having the dynamics responsible for the origin 
of the electroweak symmetry breaking separate from 
the source of the quark and lepton masses. The elec- 
troweak symmetry is assumed to be broken by some 
unspecified, presumably strong, supersymmetric dy- 
namics which further respects the custodial SU{2)v 
global symmetry. The low energy effects of the elec- 
troweak symmetry breakdown are manifested through 
an additional pair of constrained Higgs doublet super- 
fields which realize the electroweak symmetry nonlin- 
early. On the other hand, the quark and lepton su- 
perfields are Yukawa coupled only to the MSSM Higgs 
doublet superfields Hu,Hd. The additional pair of 
constrained doublet chiral superfields denoted H'^ , H'^ 



have the form H[, = 



H+' 



in+ 



with the nontrivial vacuum 

expectation values < 0\H°'\0 >^ ^ ^< 0\H^'\0 >. 
The equality of the primed VEVs is a consequence of 
the assumed SU{2)v custodial symmetry of the elec- 
troweak breaking dynamics. Implementing the super- 
field constraint, H'^eH'^ = ^, as E = ^ - H • H , 
allows the E superfield to be eliminated in favor of 
the H superfields and breaks the electroweak symme- 
try. Note that, while in the MSSM, the electroweak 
symmetry breakdown is tied to the SUSY breaking so 
that without SUSY breaking there is no electroweak 
breaking, this multi-doublet model can be realized in 
the broken electroweak symmetry phase even if SUSY 
remains unbroken. 

The Higgs doublet portion of the superpotential in- 
cludes the mixing terms among the constrained and 
MSSM Higgs muhiplets as well as the MSSM /^n-term 



so that W = fiiiHueHd + /^i2H„eff^ + fi2iH^eHd. 
The VEVs of the MSSM Higgs doublets < 0|i7°|0 >= 
< 0|i?d|0 >= are catalyzed through their bi- 
linear coupling to the constrained doublets. Note that 
nontrivial values for Vu , Vd require both nonzero v' and 
nonvanishing mixing parameters /xi2, /i2i. The quarks 
and leptons acquire masses through their couplings to 
the MSSM Higgs doublets. Finally, the soft SUSY 
breaking terms for the gauginos, A, and compo- 
nent field MSSM Higgs doublets are given by = 

\Ah (AA + AA) + \M2 {\'\' + A' A"') - mlHlHu - 
mlHlHd - fiiiBK^eHd - finBHleHl 

The model actionH] includes the SU{2)l x U{1) 
SUSY Yang-Mills action, the SUSY nonlinear sigma 
model constructed in terms of the gauged Kahler 
metric [2] -[3], the superpotential and the SUSY break- 
ing terms. After eliminating the auxiliary fields 
and imposing the constraint, the tree level poten- 
tial is constructed. There are three minimum con- 
dition which need to be satisfied. Their implemen- 
tation can be used to eliminate the SUSY break- 
ing parameters , rud and the superpotential mix- 
ing parameter /Lt2i. That leaves seven remaining 
parameters which include the MSSM mixing angle 
tan/3 — Vu/vd, the MSSM mu parameter, fin, and 
the MSSM SUSY breaking parameters Mi,M2 and 
b = —^iiB. In addition, there are two new pa- 
rameters which are the superpotential parameter ^12 
and a new mixing angle tan = (u^ -t- i'^)/2w'^. 
The total VEV is given by u = yjv^ + vl + 2v'^ 

so that Mz = \/ gl + 51^/2 = M\y/cos9w, where 
gi,g2 are the U{1), SU{2)l gauge couplings. On the 
other hand, the Yukawa coupling constants through 
which the quarks and leptons acquire their masses 
are proportionately enhanced. For example, the top 
quark Yukawa coupling is \/2mt/{v sm9 sin /3) which 
is larger than the standard model or MSSM relation. 

The stability region in parameter space is deter- 
mined by requiring all scalar squared masses to be 
positive. The model also exhibits an unbroken R- 
parity which insures the stability of the lightest SUSY 
partner (LSP). Requiring this LSP to be a neutralino 
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FIG. 1: Stability of the potential against D-flat direction runaway field values is determined in the /^ii-/ii2 parameter 
plane. Each region of SUSY breaking parameter h = —4, 000, 4, 000, 12, 000 GeV^ is depicted by the overlapping orange, 
violet, blue regions, respectively. Stability region A has tan/3 = 1, taxid = 1, region B has tan/3 = 1, tan^ = 2, region C 
has tan/3 = 2, tanS = 2, and region D has tan/3 = 10, tanS = 2. The requirement that a neutralino is the LSP further 
delineates the stability regions. The green dots indicate the points in parameters space associated with the detailed 
mass spectrum in Fig. [2] The yellow lines indicate the value of /in along which the parameter /112 is scanned in the 
subsequent mass spectrum plots. For each plot the value of the gaugino SUSY breaking masses are M\ = 200 GeV and 
M2 = 800 GeV. 
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FIG. 2: The Higgs (pseudo-) scalars and gaugino-Higgsino mass spectrum for a point in tlie LSP-stability regions 
indicated by the green dot in Fig. [l] The gaugino soft SUSY breaking masses are Mi = 200 GeV and M2 = 800 GeV, 
and b = 4,000 GeV^ for all regions. 
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FIG. 3: Masses as a function of the lightest pseudoscalar mass raa for a /ii2 scan along the yellow line across region A in 
Fig. [l] The parameters for the plots are tan/3 = 1, tan 6 = 1, 6 = 4, 000 GeV^ and /in = —12 GeV. In the top left panel 
green curves correspond to scalar h, HI, H2 masses, while the purple curve corresponds to the pseudoscalar A mass. In 
the bottom left panel, the blue curves correspond to the charged Higgs CI, C2, C3 masses. In the top right panel, the 
red curves correspond to the neutralino — masses, while the orange curves in the lower right panel correspond to 
the chargino CI, C2, C3 masses. 




FIG. 4: Masses as a function of the lightest pseudoscalar mass nria for a /ii2 scan along the yellow line across region B 
in Fig. [l] The parameters for the plots are tan/3 — 1, tanS — 2, b = 4,000 GeV'^ and /tn — —16 GeV. The curves 
correspond to the various particles just as described in the caption to Fig. [3] 
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FIG. 5: Masses as a function of the lightest pseudoscalar mass nia for a scan along the yellow line across region C 
in Fig. [l] . The parameters for the plots are tan/3 — 2, tan^ — 2, b — 4,000 GeV^ and /in — —52 GeV. The curves 
correspond to the various particles just as described in the caption to Fig. [S] 




FIG. 6: Masses as a function of the lightest pseudoscalar mass rUa for a /ii2 scan along the yellow line across region D 
in Fig. [l] The parameters for the plots are tan/3 = 10, t&nO = 2, b = 4,000 GeV^ and /in = —344 GeV. The curves 
correspond to the various particles just as described in the caption to Fig. [S] 
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further restricts the allowed parameter space. Four 
stability regions, denoted as A, B, C, and D, are 
exhibited in Fig. [T] For each panel in the figure, 
the value of the gaugino SUSY breaking masses are 
Ml = 200 GeV and Ah = 800 GeV. Stability region 
A has tan/3 = 1, tan 6* = 1, region B has tan/3 = 1, 
tan^ — 2, region C has tan/3 = 2, tan0 = 2, and 
region D has tan/3 — 10, tan 6* = 2. Each region is 
considered for three values of the SUSY breaking pa- 
rameter b = -4, 000, 4, 000, 12, 000 GeV^. In general, 
the eigenvalues of the mass matrices must be deter- 
mined numerically. Detailed mass spectra for specific 
points in parameter space indicated by green dots in 
Fig. [T] are displayed in Fig. [2j Note that the light- 
est spin zero particle can be either the neutral pseu- 
doscalar a (panels A,B) or the neutral scalar h (panels 
C, D). The next heaviest neutral pseudoscalar is de- 
noted by A, while the remaining neutral scalars in or- 
der of increasing mass are denoted as HI, H2. Adapt- 
ing a similar notation, the neutralinos in order of 
increasing mass are denoted as Nl, N2, N3, NA, N5, 
while the charged scalars (charginos) are CI, C2, C3 
(CI, C2, C3). To further explore the mass spectra, the 
neutral (pseudo-) scalar, charged scalar, neutralino, 
and chargino masses as a function of the lightest pseu- 
doscalar mass are exhibited in Figs. [3]~|6] The various 
curves in the figures follow the parameter scans from 
left to right for fixed /in with increasing /ii2 over the 
range indicated by the yellow lines in Fig. [T] for each of 
the four regions A, B, C, and D. The left endpoint of 
all the curves in each of the figures is dictated by the 
stability bounds as is the right endpoint of the curves 
in Figs. [5]|6] On the other hand, the right endpoints of 
the curves in Fig. [4] correspond to the maximum value 
for iii2 plotted in Fig. [T] All four panels allow for a 
lightest Higgs boson, h, with mass greater than 130 
GeV which is the MSSM upper bound even after the 
inclusion of radiative corrections. Using the experi- 
mental boundHj on the lightest MSSM pseudo-scalar 
of nia > 93.4 GeV as the bound for the current model, 
we see that region A allows a lightest Higgs boson tree 
level mass in the range 130 GeV < nih < 200 GeV 
which corresponds to the range 93.4 GeV < rria < 180 
GeV, while for region B, the lightest Higgs boson 
mass varies from 130 GeV < nifi < 172 GeV as 
nia ranges from 93.4 GeV < nia < 148 GeV over 
the scanned region. A lightest Higgs scalar with a 
mass in the range 115 GeV < mji < 130 GeV is 
also allowed provided different (SUSY breaking) pa- 
rameters are employed. For the scans considered, re- 
gion C admits a lightest Higgs boson mass in a range 
from 182 GeV > ruh > 115 GeV as nia varies from 
370 GeV < ma < 475 GeV. For ma less than around 
350 GeV, there is some conflict with the current ex- 
perimental limit on the mass of the lightest chargino. 
Finally region D admits a lightest Higgs boson mass 
in a range from 200 GeV > m/j > 115 GeV as ma 
varies from 3140 GeV < ma < 3180 GeV. For ma less 



than around 3000 GeV, there is some tension with the 
current experimental limit on the mass of the lightest 
chargino and/or neutralino. 

It is instructive to quantify the contributions of 
the components of the MSSM doublets and the con- 
strained Higgs doublet multiplets to the lightest de- 
grees of freedom. Here we focus on the lightest 
Higgs neutral scalar which can be decomposed as h = 
auSu + OdSd + a^S^ , where S'„ = ^{H°^ +H°),Sd = 

j^{H°J + H°) and = ^(^°t + 7^°). The vari- 
ous fractions |a„|^, |adp, |a^p are displayed in Fig. [t] 
as a function of m^. For regions A and B, a light- 
est Higgs scalar is essentially devoid of the nonlin- 
early transforming scalar St^ over the entire range 
93.4 GeV < ma- As such, the composition of the 
Higgs scalar is thus almost identical to that of the 
MSSM. In region C, the S^r fraction of is less than 
6 — 4% for a lightest Higgs scalar mass in the range 
182 GeV > mil > 115 Gev which corresponds to 
370 GeV < ma < 475 GeV. While not completely 
negligible, the Higgs scalar is still predominately com- 
posed of the MSSM fields. Finally, for region D, 
the content in the lightest Higgs scalar is about 
13 - 12% for the mass range 200 GeV > mn > 115 
GeV which corresponds to 3140 GeV < TOq < 3180 
GeV. The modification to this lightest Higgs produc- 
tion and decay due to the admixture of the non-MSSM 
content will be addressed sunsequently. The disconti- 
nuity in the slope appearing in the plots for regions A 
and _B is a consequence of the crossover in the parti- 
cle content of the lightest mass eigenvalue and the ma 
step size used in the numerical calculation. Note that 
this slope discontinuity occurs at a value of ma which 
is less than 93.4 GeV and hence excluded by the cur- 
rent experimental bound. The fractions of the lightest 
pseudoscalar, charged scalar, lightest neutralino and 
lightest chargino in terms of the MSSM fields and the 
constrained fields can be found in reference [2]. The 
field composition of the neutralino LSP is very simi- 
lar to that of the MSSM for regions A, B, C. Conse- 
quently, its identification with dark matter can pro- 
ceed just as in the MSSM. For region D, the fraction 
of TT— ino is somewhat larger being of order 10 — 5% 
for 3100 GeV < m„ < 3150 GeV. 

Finally, we briefiy address the modifications to 
Higgs boson production and decay. For moderate 
tan /3 values, the top quark loop gives the dominant 
contribution to gluon fusion Higgs production at 
the LHC provided the squark masses are sufficiently 
high 5 . The lightest Higgs boson can be writ- 
ten as a linear combination of the MSSM scalars 
Su, Sd and nonlinearly transforming scalar as 
h = UuSu + OdSd + OttStt- The modulus squares of 
various amplitudes are presented in Fig. [7] for the 
four regions of parameter space numerically probed 
in this paper. Since the top quark interacts only 
with the Su component with the enhanced Yukawa 
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FIG. 7: Lightest neutral Higgs boson, h, content as a function of the hghtest pseudoscalar mass for a ^12 scan corre- 
sponding to the yellow lines across the A, B, C, and D regions in Fig. [l] For each plot the values of the gaugino SUSY 
breaking masses are Ml = 200 GeV and M2 = 800 GeV, and 6 = 4, 000 GeV^ The scan through region A has tan /3 = 1, 
tan 6* = 1, and = —12 GeV, the one through region B has tan/3 = 1, ta.n6 — 2 and /in = —16 GeV, the one through 
region C has tan/3 = 2, tanS — 2, and /in — —52 GeV, and the one through region D has tan/3 = 10, tanO — 2, and 
/ill = —344 GeV. The red curve corresponds to the Su fraction, the green curve to the Sd fraction, and the blue curve 
to the Stt fraction. 
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FIG. 8: Ratio of gluon fusion Higgs scalar production cross-section to the standard model result. 
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FIG. 9: Ratio of two M^-boson partial decay width of the Higgs scalar, h to that of the standard model. The dashed 
line shows the enhancement (suppression) factor over the entire scanned region while solid line corresponds to the region 
where the Higgs scalar is sufficiently heavy for the decay to be kinematically allowed. 




FIG. 10: Ratio of partial width for the decay of the Higgs scalar, h, to two b quarks to that of the standard model. 



coupling -\/2m(/(usin6'sin^), the tree level gluon 
fusion production cross section is equal to that of 
the standard model times an overall factor so that 
^ = (1 + d^) (l + tisb) Note that 

the production rate depends on the details of the 
MSSM Higgs scalar Su content for the chosen values 
of parameter space. It is clear from Fig. [7] that since 
Su comprises at least one-half the Higgs scalar, there 



will be an enhanced gluon fusion production rate 
relative to the standard model as seen in Fig. [8] 
Modifications to other Higgs production processes 
such as Higgsstrahlung off a vector boson or top 
quark, or in the decay of a heavy charged Higgs 
boson, can also be considered. When considering 
the decay of the Higgs scalar, h, differences from 
the standard model arise from both the presence 
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of the two mixing angles, /?, 9, in the vacuum 
expectation values as well as the various particle 
content of h mentioned above. Since v'^ = "^'d — ^' i 
the coupling of S'tt to the W^W~ pair identically 
cancels. Consequently, the process h W^W~ 
proceeds only through the Su and Sd field compo- 
nents and the tree level decay rate of a heavy Higgs 
boson to W^W~ is the standard model rate mod- 
ified by a suppression factor so that T-^+y/- — 

(i+tan^e) (i+tan2 ,3) |a«tan/3 -I- Qd |^r™+^_ 
as seen in Fig. |9| Likewise, the decay to hb 
quarks also depends on the 6- Yukawa coupling, 
I {v sin 9 cos /3) , enhancement and the con- 
stituent fraction of the Sd content of the Higgs field. 
This leads to the modified tree level rate given by 
= |ad|2 (1 + ^) (1 + tan2 /?) T^f as displayed 
in Fig. [To] using the parameter scans appropriate 
to the four regions. For regions A and B, the 
&-pair partial rate is enhanced relative to that of the 
standard model, while for regions C and D, the rate 
is suppressed. This suppression is a consequence of 
the very small admixture of Sd in h for these regions. 

A model consisting of a supersymmetric nonlin- 
ear sigma model incorporating the low energy ef- 
fects of an unspecified electroweak symmetry break- 
ing sector and coupled to a supersymmetric version 
of the standard model was constructed and analyzed. 
The superpotential coupling of the constrained pair 
of Higgs doublets to the MSSM Higgs doublet pair 
catalyzes nontrivial vacuum expectation values in the 



later thus producing an additional contribution to the 
electroweak symmetry breaking and also serving as 
the source of the MSSM matter fields masses. The 
tree level particle spectrum of the model was obtained 
for a variety of model parameters. The MSSM upper 
limit on the mass of the lightest Higgs scalar was obvi- 
ated. Throughout the region of the explored parame- 
ter space, the lightest Higgs scalar and the neutralino 
LSP, which can be identified as a dark matter can- 
didate, were primarily composed of the MSSM fields 
with only a small admixture of the nonlinear trans- 
forming components. The main difference from the 
standard model predictions in both Higgs boson pro- 
duction from gluon fusion and Higgs scalar decay to 
either W^W^ or 66 resulted from the constituent na- 
ture of the Higgs scalar and the variant Yukawa cou- 
plings. Depending on the process and region of pa- 
rameter space, these differences could lead to either 
an enhancement or a suppression. 



Acknowledgments 

This work was supported in part by the U.S. 
Department of Energy under grant DE-FG02- 
91ER40681 (Theory and Phenomenology Task). I 
thank T.E. Clark and T. ter Veldhuis for an enjoy- 
able collaboration. 



[1] P. Fayet, Nucl. Phys. B 90 (1975) 104; J. R. El- 
lis, J. F. Gunion, H. E. Haber, L. Roszkowski and 
F. Zwirner, Phys. Rev. D 39, 844 (1989); L. Du- 
rand and J. L. Lopez, Phys. Lett. B 217, 463 (1989); 
M. Drees, Int. J. Mod. Phys. A 4, 3635 (1989); J. R. Es- 
pinosa and M. Quiros, Phys. Lett. B 279, 92 (1992); 
U. Ellwanger and M. Rausch de Traubenberg, Z. 
Phys. C 53, 521 (1992); G. L. Kane, C. F. Kolda 
and J. D. W ells, Phys. Rev. Lett. 70, 2686 (1993) 
|ar Xiv:hep-ph/9210242 ; W. T. A. ter Veldhuis, hep- 

I ph /9211281 ; J. R. Espino sa and M. Quiros, Phys. 
Lett. B 302, 51 (1993) arXiv:hep-ph/9212305 ; U. Ell- 
wanger and M. Lindner, Phys. Lett. B 301, 365 (1993) 
|arXiv:hep-ph/9211249 ; U. Ellwanger, Phys. Lett. B 
303, 271 (1993) arXiv:hep-ph/9302224 ; T. Elliott, 
S. F. King and P. L. White, Phys. Lett. B 305, 71 

(1993) arXiv:hep-ph/9302202 ; Phys. Rev. D 49, 2435 

(1994) arXiv;hep-ph/9308309 ; Phys. Lett. B 351, 213 

(1995) arXiv:hep-ph/9406303 . 

[2] T.E. Clark, S.T. Love and T. ter Veldhuis, arXiv:hep- 
ph/1107.3116. 



[3] B. Zumino, Phys. Lett. B87 (1979) 203; W. 
Buchmiiller, S.T. Love, R.D. Peccei and T. Yanagida, 
Phys. Lett. B115 (1982) 77; W. Buchmiiller, R.D. Pec- 
cei and T. Yanagida, Nucl. Phys. ^IIT (1983) 503; J. 
Bagger and E. Witten, Phys. Lett. B118 (1982) 103; 
B. Ovrut and J. Wess, Phys. Rev. D25 (1982) 409; 
T.E. Clark and S.T. Love, Nucl. Phys. B232 (1984) 
306; Phys. Lett. B 138, 289 (1984); Nucl. Phys. B254 

(1985) 569; W. Lerche, Nucl. Phys. B238 (1984) 582; 
T. Kugo, L Ojima and T. Yanagida, Phys. Lett. B135 
(1984) 402; M. Bando, T. Kuramoto, T. Maskawa and 
S. Uehara, Prog. Theor. Phys. 72 (1984) 313; 1207; K. 
Itoh, T. Kugo and H. Kunitomo, Nucl. Phys. B263 

(1986) 295; W. Buchmiiller and W. Lerche, Ann. of 
Phys. 175 (1987) 159. 

[4] K. Nakamura et. al., PDG J. Phys. G 37, 075021 
(2010). 

[5] M. Spira, A. Djouadi, D. Gra udenz and P. M. Zerwas , 
Nucl. Phys. B 453, 17 (1995) |arXiv:hep-ph/9504378] . 



